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Abstract — We investigate two quantities of interest in a delay- 
tolerant mobile ad hoc network: the network capacity region 
and the minimum energy function. The network capacity region 
is defined as the set of all input rates that the network can 
stably support considering all possible scheduling and routing 
algorithms. Given any input rate vector in this region, the 
minimum energy function establishes the minimum time average 
power required to support it. In this work, we consider a cell- 
partitioned model of a delay-tolerant mobile ad hoc network 
Vfith general Markovian mobility. This simple model incorporates 
the essential features of locaUty of wireless transmissions as 
well as node mobility and enables us to exactly compute the 
corresponding network capacity and minimum energy function. 
Further, we propose simple schemes that offer performance 
guarantees that are arbitrarily close to these bounds at the cost 
of an increased delay. 

Index Terms — delay tolerant networks, mobile ad hoc network, 
capacity region, minimum energy scheduling, queueing analysis 



I. Introduction 

Two quantities that characterize the performance limits of 
a mobile ad hoc network are the network capacity region and 
the minimum energy function. The network capacity region 
is defined as the set of all input rates that the network can 
stably support considering all possible scheduling and routing 
algorithms that conform to the given network structure. The 
minimum energy function is defined as the minimum time 
average power (summed over all users) required to stably 
support a given input rate vector in this region. Here, by 
stability we mean that the input rates are such that for all 
users, the queues do not grow to infinity and average delays 
are bounded. In this paper, we exactly compute these quantities 
for a specific model of a delay-tolerant mobile ad hoc network. 

Asymptotic bounds on the capacity of static wireless net- 
works and mobile networks are developed by |2|, |3|. The 
work in p| shows that for networks with full uniform mobility, 
if delay constraints are relaxed, a simple 2-hop relay algorithm 
can support throughput that does not vanish as the number of 
network nodes grows large. Recent work in |]4) generalizes 
this model and investigates capacity scaling with non-uniform 
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node mobility and heterogeneous nodes. Capacity-delay trade- 
offs in mobile ad hoc networks are considered in ISl-fl^. 
Flow-based characterization of the network capacity region is 
presented in several works (e.g., ||7), p3) , p4)). 

However, little work has been done in computing the 
exact capacity and energy expressions for these networks. 
Exceptions include a closed form expression for the capacity 
of a fixed grid network in |5|, an expression for the exact 
information theoretic capacity for a single source multicast 
setting in a wireless erasure network ||6), and an expression for 
the capacity of a mobile ad hoc network in |8 | that uses a cell- 
partitioned structure. The work in |8| quantizes the network 
geography into a finite number of cells over which users move, 
and assumes that a single packet can be transmitted between 
users who are currently in the same cell, while no transmission 
is possible between users currently in different cells 

In this work, we extend this model to more general scenar- 
ios allowing adjacent cell communication and different rate- 
power combinations. Specifically, we extend the simplified 
cell-partitioned model of |8| (which only allows same cell 
communication and considers i.i.d. mobility) to treat adjacent 
cell communication. We establish exact capacity expressions 
for general Markovian user mobility processes (possibly non- 
uniform), assuming only a well-defined steady-state location 
distribution for the users. Our analysis shows that, similar 
to |8J, the capacity is only a function of the steady-state 
location distribution of the nodes and a 2-hop relay algorithm 
is throughput optimal for this extended model as well. Further, 
our analysis illuminates the optimal decision strategies and 
precisely defines the throughput optimal control law for choos- 
ing between same cell and adjacent cell communication. We 
then use this insight to design a simple 2-hop relay algorithm 
that can stabilize the network for all input rates within the 
network capacity region. We also compute an upper bound on 



the average delay under this algorithm. (Sec. Ill i 

We next compute the exact expression for the minimum 
energy required to stabilize this network, for all input rates 
within capacity. Our result demonstrates a piecewise linear 
structure for the minimum energy function that corresponds 
to opportunistically using up successive transmission modes. 
Then we present a greedy algorithm whose average energy 
can be pushed arbitrarily close to the minimum energy at the 
cost of an increased delay. (Sec. |1V|| 



' Here a cell represents only a sub-region of the network. There are no base 
stations in the cells (i.e., this should not be confused with "cellular networks" 
that use base stations.) 
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Fig. 1 . An illustration of the cell-partitioned network with same and adjacent 
cell communication. Cells that share an edge are assumed to be adjacent. 

Before proceeding further, we emphasize that the network 
capacity and minimum energy function derived in this paper 
are subject to the scheduHng and routing constraints of our 
model as described in the next section. Specifically, in this 
work, we do not consider techniques that "mix" packets, such 
as network coding or cooperative communication, which can 
increase the network capacity and reduce energy costs. In fact, 
in Sec. [V] we present an example scenario that shows how 
network coding in conjunction with the wireless broadcast 
advantage can increase the capacity for this model. Calculating 
the network capacity region and the minimum energy function 
when these strategies are allowed is an open problem in 
general in network information theory and is beyond the scope 
of this paper. 

II. Network Model 

A. Cell-Partitioned Structure 

We use a cell-partitioned network model (Fig. [TJ having C 
non-overlapping cells (not necessarily of the same size/shape). 
There are N users roaming from cell to cell over the network 
according to a mobility process. Each cell c G {1, 2, . . . . C} 
has a set of adjacent cells Be that a user can move into from 
cell c. The maximum number of adjacent cells of any cell is 
bounded by a finite constant J. We define the network user 
density ■&& d ~ N/C users/cell. For simplicity, N is assumed 
to be even and N > 2. Note that there could be "gaps" in 
the cell structure due to infeasible geographic locations. We 
assume that the gaps do not partition the network, so that it 
is possible for a single user to visit all cells. We assume C is 
the number of valid cells, not including such gaps. 

B. Mobility Model 

Time is slotted so that each user remains in its current 
cell for a timeslot and potentially moves to an adjacent cell 
at the end of the slot. We assume that each user i moves 
independently of the other users according to a mobility 
process that is described by a finite state ergodic Markov 
Chain. In particular, let P = {Pij}cxc be the transition 
probability matrix of this Markov Chain. Then Pij represents 
the conditional probability that a user moves to cell j in the 
current slot given that it was in cell i in the last slot. Note 



that Pij > only if j is an adjacent cell of i, i.e., j e Bj. It 
can be shown that the resulting mobility process has a well- 
defined steady-state location distribution tt = {tTc}ixc over 
the cells c e {1,2,..., C} that satisfies ttP = tt and is the 
same for all users. However, this distribution could be non- 
uniform over the cells. We assume that in each slot, users are 
aware of the set of other users in the same cell and in adjacent 
cells. However, the transition probabilities associated with the 
Markov Chain P are not necessarily known. 

It can be shown (see, for example, p8) ) that the mobility 
process discussed above has the following property. Let x(i) e 
{1, . . . , C} denote the location of a user in timeslot t. Then, 
for all integers d > 0, there exist positive constants a, 7 such 
that Vc e {1, 2, . . . , C}, the following holds: 

TTcil - ai"") < Pr[x{t + d)^ c\x{t)] < n,{l + aj"") (1) 

where a > 1 and < 7 < 1. Moreover, the decay factor 
7 is given by the second largest eigenvalue of the transition 
probability matrix P (see |17|). From this, it can be seen 
that for any e > 0, choosing d = [tI^^I ensures that the 
conditional probability that the user is in cell c at time t + d 
is within TTce of the steady-state probability tTc of being in 
cell c, irrespective of the current location. This implies that 
the Markov Chain converges to its steady-state probability 
distribution exponentially fast. Using the independence of 
user mobility processes, the following can be shown about 
functionals of the joint user location process xi^)- 

Lemma 1: Let x{t) — {xi{t)i ■ ■ ■ iXN{t)) be the vector of 
current user locations, where Xi(i) represents the cell of user 
i in slot t. Let f{x{t)) be any non-negative function of x(^)> 
i.e., f{x{t)) > V x{t)- Define fav as the expectation of 
/(x(t)) over the steady-state distribution of x{t)- 

N 

fav^ X! /(ci,---,Cjv) JItTc, 

Ci,C2,...,Cjv i=l 

Then for all d such that aj"^ < j^j, we have 

fa.il - 2Naj'') <E{f{x{t + d))\x{t)} < /a.(l + 2Naj'i) 

Proof: See Appendix A. ■ 

C. Traffic Model 

We assume that there are N unicast sessions in the network 
with each node being the source of one session and the 
destination of another session. Packets are assumed to arrive 
at the source of each session i according to an i.i.d. arrival 
process Ai{t) of rate A^. We assume that in any slot, the 
maximum number of arrivals to any session i is bounded, 
i.e., Ai{t) < Ajnax- While our analysis holds for the general 
source-destination pairing, for simplicity, we assume that N 
is even with the following one-to-one pairing between users: 
1 O 2, 3 O 4, . . . , (A^ - 1) O A^, i.e., packets generated by 
user 1 are destined for user 2 and those generated by user 2 
are destined for user 1 and so on. This assumption simplifies 
the computation of the capacity in closed form in Theorem [T| 
and will be used for the rest of the paper 
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D. Communication Model 



III. Network Capacity 



We assume that two users can communicate only if they 
are in the same cell or in adjacent cells. Further, if the 
communication takes place in the same cell, i?i packets can be 
transmitted from the sender to the receiver if the sender uses 
full power If the receiver is in an adjacent cell, R2 packets 
can be transmitted with full power We assume that Ri and R2 
are non-negative integers and that Ri > i?2. Power allocation 
is restricted to the set {0, 1}, i.e., each user either uses 
zero power or full power For simplicity, we assume that the 
communication cost consists only of the transmission power. 
The analysis presented can be easily extended to the case 
with non-zero reception power by defining the communication 
cost as the total power (including transmission and reception) 
required for sending -Ri(i?2) packets from a transmitter to a 
receiver in the same (adjacent) cell. 

We allow at most one transmitter in a cell at any given 
time slot, though the cell may have multiple receivers (due 
to possible adjacent cell communication). Further, a user may 
potentially transmit and receive simultaneously. This model 
is conceivable if the users in neighboring cells use orthogo- 
nal communication channels. This model allows us to treat 
scheduling decisions in each cell independently of all other 
cells, thereby enabling us to derive closed form expressions 
for capacity and minimum energy. 



E. Discussion of Model 

While an idealization, the cell-partitioned model captures 
the essential features of locality of wireless transmissions 
as well as node mobility and allows us to compute exact 
expressions for the network capacity and minimum energy 
function. This model is reasonable when nodes use non- 
interfering orthogonal channels in adjacent cells. We also 
refer to Section I-A of jSJ for further discussion on the cell- 
partitioned network assumption. 

In this work, we restrict our attention to network control 
algorithms that operate according to the network structure 
described above. A general algorithm within this class will 
make scheduling decisions about what packet to transmit, 
when, and to whom. For example, it may decide to transmit to 
a user in an adjacent cell rather than to some user in the same 
cell, even though the transmission rate is smaller. However, 
we assume that the packets themselves are kept intact and are 
not "mixed" (for example, using cooperative communication 
and/or network coding). Allowing such strategies can increase 
the capacity, although computing the exact capacity region 
remains a challenging open problem in general. In Sec. |V] 
we present an example that shows how network coding in 
conjunction with the wireless broadcast advantage can increase 
the capacity for this model. However, we note that if we 
remove the broadcast feature, then the scenario considered in 
this paper becomes a network coding problem for multiple 
unicasts over an undirected graph, for which it is not yet 
known if network coding provides any gains over plain routing 
(see further discussion in |16|). 



In this section, we compute the exact capacity of the net- 
work model described in the previous section. For simplicity, 
we assume that all users receive packets at the same rate (i.e., 
Xi = A for all i). The capacity of the network is then described 
by a scalar quantity which denotes the maximum rate A that the 
network can stably support. Recall that network user density 
9 = N/C users/cell. Then we have the following: 

Theorem 1: The capacity of the network (in packets/slot) is 
given by: 



{ Riq+B.ip+R2q'+B.2p' 
2'?, „ , , 

2Riq+2R2q"+Rip"+B.2{p -q 
26 



if i?i > 2R2 

if 2i?2 > Ri>R2 



^ X^^sLi i^J'ifinding no source-destination pair in cell c 



where 

q — ^Y^c=i ^''^ [finding a source-destination pair in cell c in 
a timeslot] 

P — Y^c=i ^'J' [finding at least 2 users in cell c in a timeslot] 
q' — ^ X^^jLi [finding exactly 1 user in cell c and its 
destination in an adjacent cell in a timeslot] 
p' = ^ X]c=i J^'' [finding exactly 1 user in cell c and at least 
1 user in an adjacent cell in a timeslot] 

<l"-c " 

but at least 1 source-destination pair with an adjacent cell in 
a timeslot] 

p" = ^ Y^c=i ^'''[finding no source-destination pair in cell c 
and any adjacent cell but at least 2 users in the cell c in a 
timeslot] 

The probabilities in the summations above are the proba- 
biUties associated with the steady-state cell location distribu- 
tions of the users. Using the independence of user mobility 
processes and the same steady-state cell location distribution 
TT = {tTcIixc for all users, we can exactly compute these 
probabilities for our model. These are given by (see Appendix 
B for detailed derivation): 



c=l 
C 



P = 



c— 1 



(1 - (1 - vr,)^ - Ntt,{1 - vr,)^-!) 
^^(n„rf,(c)iV7r,(l-7r,)^ 



\N-1 



c=l 
C 



n ^ 

c=l i=l 
n ^ 

."-^EE2^ 



(1 - (1 - n,,, (c))~-i)7V7r,(l - ttJ^-i 



c=l 



7r:(i-^,)^-Mi-(i-n,,,(c))' 



(2) 



Here, Iladj{c) denotes the sum of the conditional steady-state 
probability of a user being in any adjacent cell of cell c given 
that this user is not in cell c, i.e., Iladj{c) = i\ Tli^B "^i- 
Thus, the network can stably support users simultaneously 
communicating at any rate A < /i. We prove the theorem 



4 



in two parts. First, we establish the necessary condition by 
deriving an upper bound on the capacity of any stabilizing 
algorithm. Then, we establish sufficiency by presenting a 
specific scheduling strategy and showing that the average delay 
is bounded under that strategy. 



A. Proof of Necessity 

Proof Let "if be the set of all stabilizing scheduling 
policies. Consider any particular policy i/; e Vf. Suppose it suc- 
cessfully delivers (T) packets from sources to destinations 
involving "a" same cell transmissions and "&" adjacent cell 
transmissions in the interval (0,T). Fix e > 0. For stability, 
there must exist arbitrarily large values of T such that the total 
output rate is within e of total input rate. Thus: 



of packet transmissions (for example, when enough packets are 
not available). 



Let Z-(r) - r-(r) + l^(,^,(r) + X^.Jt). Also define 
the following indicator decision variables for any policy ui for 
some T e (0, T) and c e {1, 2, . . . , C}: 



ECO \— ^ x: v'V' 



(T) 



T 

as the total 



> iVA 



(3) 



Define Y'^{T) as the total number of packet transmis- 
sions in (0,T) under pohcy ip. Then, V^iT) is at least 
J2T=o '^'bLoi'^~^^)-^tbi'^) (because these many packets were 
certainly delivered). Thus, we have 



1 

f' 



Y^iT)>-J2J2i^ + b)Xt,iT) 



> 



> 



a=0 fc=0 



E 

a+b<2 

E 

a+b<2 



XtbiT) + 



E 

a+b>2 



XtbiT) 



Xt,{T) + 2{N\-e)- 



T 



E ^afc 

-i+b<2 



(T) 



where the last inequality is obtained using ([3]l. Hence, noting 
that e can be chosen to be arbitrarily small, we have: 



A < lim 

T-i-oo 



xUt) 



^01 



(T) 



2TN 



(4) 



Define Y^{t) as the total number of packet transmissions in 
cell c at timeslot r under policy ip. Also define xfQ ^(r) and 
-Xqi c(^) number of packets delivered by same cell 



direct and adjacent cell direct transmission respectively in cell 
c at timeslot t. Then y'^(T) + xfg{T) + X^^{T) can be 
written as a sum over all timeslots t e (0, T) and all cells c 
as follows: 



Y^iT)+xUT)+X;},iT) 
= E E (Yc^'ir) + Xtojr) + Xtar) 



T = C=l 

T-1 C 

^ E E "^ii i^cir) + X-o,,{r) + X-,^,{r)) (5) 

where Y^ (t) denotes the total number of packet transmission 
opportunities in cell c at timeslot t under any policy u. 
Similarly, Xi^ ^(r) and X^i ^{t) denote the total number of 
packet transmission opportunities that use same cell direct 
and adjacent cell direct transmissions respectively in cell c 
at timeslot t. Note that these do not depend on the queue 
backlogs and therefore can be different from the actual number 



m = 



Itir) = 



1 if a same cell direct transmission can 
be scheduled in cell c in slot r 

else 

1 if a same cell relay transmission can 
be scheduled in cell c in slot r 

else 

1 if an adjacent cell direct transmission can 
be scheduled in cell c in slot r 

else 

1 if an adjacent cell relay transmission can 
be scheduled in cell c in slot t 

else 



Note that the transmission rates associated with these decision 
variables are i?i , i?i , R2 and R2 respectively. Then, we can 
express Z'^{t) as follows: 

i-(r) = 1;"(t) + Xro,,(r) + X^.Ar) = Ril'dr) + Rill{r) 
+ R2ll{T) + R^ltir) + X- ,,(r) + ,,(r) 
= RiII{t) + RiII{t) + i?2/c (r) + i?2/c (r) + RiII{t) 
+ i?2/c (r) = 2i?i/i(r) + RiII{t) + 2i?24^(r) + i?2/c (r) 

Note that under any scheduling policy, only one of the decision 
variables can be 1 and the rest are 0. Thus, the preference order 
for decisions to maximize Z'^{t) is evident. Specifically, it 
would be /1(t) )~ I^{t) >- I fir) /^(t) when i?i > 2i?2 
and /1(t) y I^{t) >- I^{t) y I^{t) when R2 < Ri < 2i?2. 
Thus, in each cell c, Z^{t) is maximized by the policy oj 
that chooses the scheduling decisions in this preference order, 
choosing a less preferred decision only when none of the more 
preferred decisions are possible in that cell. 



Define Zc{t) = max^g,j, Z"(t). Then using Q and (U, 
we have 

T-l C 

A< lim T^T^rYY Z,{t) 

T = C=l 

As Zc(r) can take only a finite number of values (namely 
i?i , i?2 , 2i?i , 2i?2 and 0) and is a function of the current state 
of the ergodic user location processes, the time average of 
Zc(t) is exactly equal to its expectation with respect to the 
steady-state user location distribution. Thus, the bound above 
can be computed by calculating the expectation of Z(.(t) using 
the steady-state probabilities associated with the indicator 
variables and summing over all cells. When R\ > 2i?2, this 
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is given by: 



1 



T-l C 

C 



^EE{^c(r)} 



2N 



c=l 



2Riq + Ri{p-q) + 2R2q' + R2{p' - q') 
20 



and when R2 < Ri < 2i?2, this is given by: 
1 



T-l C 
r=0 c=l 



^EE{Z,(r)} 



2iV 

c=i 

^ 2i?ig + 2^2 + Rip" + i?2(y - q') 
26 

This establishes the necessary condition for the network 
capacity. ■ 

Note that the above preference order clearly spells out 
the structure of the throughput optimal strategy. Specifically, 
depending on the values of Ri and R2, this order can be 
used to decide between same cell relay and adjacent cell 
direct transmission. We use this insight to design a throughput- 
optimal 2-hop relay algorithm in the next section. Also note 
the factor of 2 with the decision variables corresponding to 
direct source-destination transmission. Intuitively, each such 
transmission opportunity is better than a similar opportunity 
between source-relay or relay-destination by a factor of 2 since 
the indirect transmissions need twice as many opportunities 
to deliver a given number of packets to the destination as 
compared to direct transmissions. 



B. Proof of Sufficiency 

Now we present an algorithm that makes stationary, ran- 
domized scheduling decisions independent of the queue back- 
logs and show that it gives bounded delay for any rate A < /i, 
i.e., there exists a p such that < p < 1 and A — pfi. We only 
consider the case when i?i > 2i?2. The other case is similar 
and is not discussed. 

2-Hop Relay Algorithm: Every timeslot, for all cells, do 
the following: 

1) If there exists a source-destination pair in the cell, ran- 
domly choose such a pair (uniformly over all such pairs in 
the cell). If the source has new packets for the destination, 
transmit at rate Ri. Else remain idle. 

2) If there is no source-destination pair in the cell but there 
are at least 2 users in the cell, randomly designate one 
user as the sender and another as the receiver. Then, with 
probability (where {) < 5 < 1 and S = 6{p) to be 
determined later), perform the first action below. Else, 
perform the second. 

a) Send new Relay packets in same cell: If the transmitter 
has new packets for its destination, transmit at rate 
Else remain idle. 



b) Send Relay packets to their Destination in same cell: 
If the transmitter has packets for the receiver, transmit 
at rate Else remain idle. 

3) If there is only 1 user in the cell and its destination is 
present in one of the adjacent cells, transmit at rate R2 
if new packets present. Else remain idle. 

4) If there is only 1 user in the cell and its destination is not 
present in one of the adjacent cells but there is at least 
one user in an adjacent cell, randomly designate one such 
user as the receiver and the only user in the cell as the 
transmitter. Then, with probability (where < 5 < 1 
and (5 = 5{p) to be determined later), perform the first 
action below. Else, perform the second. 

a) Send new Relay packets in adjacent cell: If the trans- 
mitter has new packets for its destination, transmit at 
rate i?2. Else remain idle. 

b) Send Relay packets to their Destination in adjacent 
cell: If the transmitter has packets for the receiver, 
transmit at rate i?2. Else remain idle. 

This algorithm is motivated by the proof of necessity of 
Theorem [T] since it follows the same preference order in 
making scheduling decisions. Note that this algorithm restricts 
the path lengths of all packets to at most 2 hops because any 
packet that has been transmitted to a relay node is restricted 
from being transmitted to any other node except its destination. 

To analyze the performance of this algorithm, we make 
use of a Lyapunov drift analysis |7|. Consider a network 
of N queues operating in slotted time, and let U{t) = 
{Ui{t), U2{t), . . . , C/Ar(t)) represent the vector of backlogged 
packets in each of the queues at timeslot t. Let L{U{t)) be 
a non-negative function of the unfinished work U{t), called 
a Lyapunov function. Define the conditional Lyapunov drift 
A{t, d) at time t > d (where d > in a fixed integer) as 
follows: 

A{t, d)4E {L{U{t + 1)) - L(U{t))\U{t - d)} 

Then we have the following lemma. 

Lemma 2: Lyapunov Drift Lemma: If there exists a positive 
integer d such that for all timeslots t > d and for all U {t), the 
conditional Lyapunov drift A{t,d) satisfies: 



N 



A{t,d) <B-e^U,{t~d) 
1=1 



(6) 



for some positive constants B and e, and if E |L([7(d))| < 
00, then the network is stable, and we have the following 
bound on the time average total queue backlog: 

t-l N 



limsup^ VVE{C/,(r)} 



< 



B 



(7) 



r=0 1=1 



Proof This can be shown using a telescoping sum 
technique (similar to related proof in |7|) and is omitted for 
brevity. ■ 

We now make use of this lemma to analyze the performance 
of the 2-Hop Relay Algorithm. 

Theorem 2: For the cell partitioned network (with N nodes 
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and C cells) as described in Sec. |ll] with capacity fi = 

Rip+Rip'+Riq+R^q' jjjj^J jjjp^j jj^jg )^ g^gjj ^ggj. g^gjj jjjj^j 

A = p/i for som e < p < 1, and user mobility model as 
described in Sec. II-B the average packet delay D under the 
2-Hop Relay Algorithm with 5 — ^-r^ satisfies: 



D < 



BN{2d+l) 



(8) 



XflK{l — p) 

where B is a constant given by ( [TT] ), k is a positive constant 

B'''tB'^To"^7B"', , and d is a finite integer 

R1P+R2P' +Riq+R2q ^ 



given by k 

that is related to the mixing time of the joint user mobility 

- i°g(^) " 

log(l/7) 



process and is given by d = 



Proof: Let u\'^\t) represent the total backlog of type c 
(i.e., number of packets destined for node c) that are queued 
up in node i at time t. The queueing dynamics of ul'^\t) 
satisfies the following for all c ^ i: 



C/f^(i + 1) < max Ul^'it) 



b 



(9) 



where A'f' {t) =number of new type c arrivals to source node 
i at the beginning of timeslot t and ^^l^^{t) =rate offered to 
type c packets in timeslot t with node a as transmitter and 
node b as receiver. The above is an inequality because the 
actual number of packets transmitted from the other nodes 
to node i (for relaying) could be less than the incoming 
transmission rate p!^^}^ (t) when these nodes do not have 
enough packets. Now define a Lyapunov function L{U{t)) — 
T!i'=iY.c^i{Ui''\t)Y ■ Using the conditional Lyapunov 
drift A(t, d) can be expressed as follows: 



N 



= 1 c^i 



(: 



b a 

Here, B is given by: 

B = {A 



{t)~A^\t))\U{t-d) 



out \ 
max ) 



(10) 



(11) 



where p^ax — maximum transmission rate into any node = 
i?i + JR2, where J is the maximum number of adjacent cells 
of any cell (Sec. II-Ai and — maximum transmission 



rate out of any node ~ Ri. 

We now use the following sample path relations to express 
(10 1 in terms of the queue backlog values at time {t — d). 



Specifically, we have the following for alH > d where d is a 
positive integer (to be determined later) for all i ^ c. 

J2 Ut\t -d)+ d{A^ax + fCa.) > E 

Y^ui^\t-d)-dpi:i<Y^ui^\t) 

These follow by noting that the queue backlog at time t cannot 
be smaller than the queue backlog at time [t — d) minus the 
maximum possible departures in duration {t — d, d). Similarly, 



it cannot be larger than the queue backlog at time {t — d) plus 
the maximum possible arrivals in duration {t — d,d). Using 
these, we can express ( [TO) i in terms of the "delayed" queue 
backlogs f/j^^' {t — d) as follows: 

N 

A{t, d) < BN{2d + 1) - 2 ^ ^ C/f ^ (t - d) X 



i — 1 C^i 



E 



Af>{t)\U{t^d) 



(12) 



Let T(t — d) = {x{t — d),U{t — d)) represent the composite 
system state at time (t — d) given by the user locations and 
queue backlogs. Since the 2-Hop Relay Algorithm makes con- 
trol decisions only as a function of the current user locations, 
the resulting service rates are functionals of the Markovian 
mobility processes. By the Markovian property of the x{t~d,) 
process, any functionals of this also converge exponentially 
fast to their steady-state values. Thus, using Lemma [T] when 
aj"^ < l/N'^, we have the following bounds: 



E 



Y.P^\t)\U(t-d) 



f^'-,:^it)\Uit-d) 



> 



< 



E 

6 

E 



JI^^)il-2Na^'') (13) 
pif)(l + 2iVa7'') (14) 



where fJ-li, , fJ-g^/ are the steady-state service rates achieved by 
the 2-Hop Relay Algorithm. We now compute these values 



and use the inequalities (13 1, ( 14 1 to obtain a bound on ( 12 1. 
We have the following 2 cases: 

1) Node i Generates Type c Packets: In this case, 
E|A|'=^(t)| = A and EaMi^V) = (since under the 2- 
Hop Relay Algorithm, a source node would never get back a 
packet that it generates). To calculate X^hMib^ we note that 
the outgoing service rate for packets generated by the source 
is equal to the sum of the rate at which the source is scheduled 
to transmit directly to its destination and the rate at which it is 
scheduled to transmit type c packets to any of the relay nodes. 
Let these rates be ri and r2 respectively. Also let the trans- 
mission rate at which it is scheduled to transmit relay packets 
to their destinations be r-^. Since the 2-Hop Relay Algorithm 
only schedules transmissions of these types, the total rate of 
transmissions over the network is given by N{ri + r2 + ra). 
Using the probability of choosing source-relay and relay- 
destination transmissions, we have: r2 = ifpf'^s- In the 2-Hop 
Relay Algorithm, a direct source-to-destination transmission 
is scheduled whenever there is a source-destination pair in the 
same cell or there is only 1 node in a cell and its destination 
is in an adjacent cell (and independent of the actual queue 
backlog values). Thus, using the definitions of q and q' from 
the statement of Theorem[T| we have: Nri — C{Riq + R2q')- 
Similarly, the sum total transmissions in the network can be 
expressed in terms of the quantities p and p' as follows: 
N{ri + r2 + 7-3) — C{Rip + i?2P')- Using these to solve 
for ri, r2, r3 and simplifying, we have 

ri=/^(l-K), r2 = ^i^^{l-5), = ^k{1 + 5) (15) 
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where n=^r^^f^TT^^^^^ ■ Note that < k < 1 (since 
p > q and p' > q'). Therefore, we have: 



b 



ri + r2 = /I - plk6 



Let 6 = ^ and a^'^ = 



5 _ i-p 

(5 represents a vahd probability since < p < 1 



Note that this choice of 
< p < 1. Then, using 
( [T3] l, the last term of ( [T2| under this case can be expressed as: 



b 



where we used the fact that (1 - kS){1 - ^) > (1 - S)^. 

2) Node i Relays Type c Packets: Note that TV > 2 for 
this case to happen. From our traffic model, we know that 
in this case A'f\t) — for all t. Further, under the 2-Hop 
Relay Algorithm, /i^^' (t) > only if node a is the source for 
type c packets. Also /i^Jj' (t) > only if 6 c. To compute 

SbMib'' ^nd J^aT^aA '^his case, note that the 2-Hop Relay 
Algorithm schedules relay transmissions such that all {N — 2) 
relay packet types are equally likely. Thus we have: 

rs -(c) r2 



(c) 
ib 



3 

—0^ E^ 



iV-2 



iV- 2 



Let 5 



V and a7'^ = 2^ = 1^. Then, using (Tl, (14), 

isea as: 



the last term of (12 1 under this case can be expressed 



E I E ^^ (*) - E ^^ (^) - ^ (^) \u{t - d) I 

> ( E T^^Y^ - 27Va7'^) - ( E Pif ) (1 + ^Na^'^) 



b a b a 

(r3 - r^) - ^^^^ ^ ^ / _ - P) 

A^-2 N -2\ n) - 2N 



where we used ilSi. 

and a-f'^ = ^0: 



Combining these two cases, with 6 = -^-j^ 



E \T.^^(t) -E A^if w -^'^(t - d)] > 



Using this in (12i, we get 



A{t,d) < BN{2d+l) - 



N 



N 



EY.ur\t-d) 

2 — 1 C^i 



This is in a form that fits (|6|. Using the Lyapunov Drift 
Lemma, we get 



BN^{2d+l) 



T—0 i^c 



Hk{1 - p) 



(16) 



^x = 0.05 
— x = 0.1 
— x = 0.25 

analytical bound 
for l.l.d. mobility 




0.2 0.25 0.3 0.35 
input rate X (packets/slot) 



Fig. 2. Average packet delay under the 2-Hop Relay Algorithm in a network 
of 16 cells with 20 nodes for different mixing times of the mobility process. 



The total input rate into the network is NX. Thus, using 
Little's Theorem, the average delay per packet is bounded by 

BN{2d+l) ^ 

C. Discussion and Simulation Example 

The proof of the capacity for the cell-partitioned network 
can be used to consider more general scheduling restrictions. 
From it amounts to: 



1 



c 



A < — — E < max / 

\ c=l 



If the bound on the right hand side can be achieved by any 
policy (potentially randomized) that takes decisions only as a 
function of the current network state, then we can design a 
deterministic policy that is throughput optimal by scheduling 

to maximize Y!I=i {Yc{i) + X^oA^) + ^oi,c(^)) subject to 



the network restrictions. For the specific cell-partitioned model 
considered here, this maximization is achieved by following 
the preference order of the decision variables in each cell 
separately as described earlier. This enables us to exactly 
compute the capacity of the network. It is possible to do the 
same for extensions to this model involving other constraints. 
For example, under the constraint that a user cannot simultane- 
ously transmit and receive, the above maximization becomes a 
maximum-weight match problem. Similarly, one could allow 
more than one transmitter per cell, in which case we would 
need to define more indicator decision variables for all possible 
control options. 

We next consider an example network consisting of 20 
nodes and 16 cells as shown in Fig. [T] The nodes move from 
one cell to another independently according to a Markovian 
random walk. Specifically, at the end of every slot, a node 
stays in its current cell with probability (1 — x), else it 
decides to move randomly one step in either the North, West, 
South, or East directions with probability x. If there is no 
feasible adjacent cell, then the user remains in the current 
cell. It can be shown that the resulting steady-state location 
distribution is uniform over all cells for all < x < 1. Thus, 
T^c = yq for all cells c. Next we assume that Ri = 2 and 
i?2 = 1 packets/slot. Then using Theorem [T] the capacity for 
this network is given by = ^^^i+Rip+^i'+R^p and can 
be calculated exactly. Specifically, we get p — 0.358, q — 
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0.038, p' = 0.357, q' — 0.073 and the network capacity is 
given hy fi — 0.489 packets/slot. 

We next simulate the 2-Hop Relay Algorithm on this 
network. New packets arrive at each source node according 
to independent Bernoulli processes, so that a single packet 
arrives i.i.d. with probability A every slot. In Fig. |2j we plot 
the average packet delay vs. A for different values of x. We 
also plot the analytical bound (j8]l of Theorem |2] for the i.i.d. 
mobility case (for which d = 0). It can be seen that the 
average delay goes to infinity as A is pushed closer to the 
capacity = 0.489 packets/slot (shown by the vertical line 
in Fig. |2]l. While the network capacity is the same for all 
values of x (since x does not affect the steady-state location 
distribution), the average delay increases as x becomes smaller 
This is because a smaller x implies a larger value for the 
parameter d leading to larger delay as suggested by the delay 
bound ([8]) in Theorem |2] Thus, the 2-Hop Relay Algorithm 
is able to support all input rates within the network capacity 
with finite average delay. However, its delay performance is 
not necessarily the best. For example, when the input rate is 
small (say A = 0.1 packets/slot), the average delay is more 
than 100 slots. Note that the 2-Hop Relay Algorithm makes 
scheduling decisions purely based on the current user locations 
and restricts all packets to at most 2 hops. It does not attempt 
to optimize the delay in the network. The delay performance 
may be improved using alternative scheduling strategies that 
do not restrict packets to at most 2 hops. For example, backlog 
aware scheduling and routing (e.g., or schemes that exploit 
the mobility pattern of the users (e.g., ifTSl) may offer better 
delay performance. 

IV. Minimum Energy Function 

We now investigate the minimum energy function of the 
cell-partitioned network under consideration. Recall that in our 
network model, each user either uses zero power or full power 
Further, i?i(i?2) packets can be transmitted from the sender 
to the receiver in the same (adjacent) cell if the sender uses 
full power 

The minimum energy function <1>(A) is defined as the 
minimum time average energy required to stabilize an input 
rate A per user, considering all possible scheduling and routing 
algorithms that conform to the given network structure. We 
exactly compute this function for our network model. Specif- 
ically, we assume that all users receive packets at the same 
rate (i.e., Ai — A for all i). Also, we consider the case when 
Ri > 2i?2 (<i>(A) for the case when i?i < 2i?2 has a different 
expression, but the proof is similar). 

Theorem 3: The minimum energy function i'(A) per user 
for the cell-partitioned network as described in Sec. |ll] with 
Ri > 2i?2 is a piecewise linear curve given by the following: 



Riq 



— Riq 



< A < 



Rijp+q) 



C. EE 



$(A) = < 




2-B.19 
29 

Rl {P+<1) \ 
26 J 

Rl(p+q)+2R2q' 
26 



A- 



if Cl 
if C2 
if C3 

if C4 



where Ci = 0<X<^,C2=-^:^a-^ - - , 03 

Rl{P+q) / \ ^ Ri(p+q)+2R2q' ^ _ Ri(p+q)+2R2q' / \ ^ 
26 ^ < 20 '<-"4 = 26 S A < 

/i. Thus, the network can stably support users simultaneously 
communicating at any rate X < fi with an energy cost that can 
be pushed arbitrarily close to $(A) (at the cost of increased 
delay). We prove the theorem in two parts. First, we establish 
the necessary condition by deriving a lower bound on the 
energy cost of any stabilizing algorithm. Then, we establish 
sufficiency by presenting a specific scheduling policy and 
showing that the average delay is bounded under that policy. 

A. Proof of Necessity 

Proof: Consider any scheduling strategy that stabilizes the 
system. Let Xab{T) denote the number of packets delivered 
by the strategy from sources to destinations in time interval 
(0,T) that involves exactly a same cell and b adjacent cell 
transmissions. For simplicity, assume that the strategy is er- 
godic and yields well defined time average energy expenditure 
per user e and well defined time average values for Xab where: 

Xab{T) 



lim 



The average energy cost per user e of this policy satisfies: 

b \ Xab 

R2) N 



(17) 



e > 



y(- 

a.b 



(18) 



This follows by noting that enough packets may not be 
available during a transmission. 

Note that xqo = 0, and so the only possible non-zero Xab 
variables are for (a, b) pairs that are integers, non-negative, 
and such that (a, 6) 7^ (0,0). Let x — (xab) represent the 
collection of Xab variables, and note that these variables must 
satisfy the constraint a; G fio H fii n il2 H fl^, where the four 
constraint sets are defined below: 



n2^< 



(a,&)#(0,0) 



Xab = NX 



0i4<^ X 



xio 
Rl 



< Cl 



1 

R~i 



aXaO < Cl + C2 



-E 



aXaO 



R2 



< Cl + C2 + C3 



where ci is the maximum rate of source-destination trans- 
mission opportunities in the same cell, ci + C2 is the maxi- 
mum rate of all possible same cell transmission opportunities 
and Cl + C2 + C3 is the maximum rate of all same cell 
or source-destination adjacent cell transmission opportunities. 
Here, these quantities are summed over all cells. Using the 
definitions of p, q and q' from the statement of Theorem [T] we 
know that Cl = Cq,ci + C2 = Cp, C1+C2+C3 — C{p+q'). For 
example, (ci + C2 + C3) can be written as ^ J2l=o (^i(^) + 

X2{t)+X3{t)j where Xi{t) is the maximum number of direct 
same cell opportunities, X2{t) is the maximum number of 
indirect same cell opportunities given all direct opportunities 
are used and X3(t) is the maximum number of direct adjacent 
cell opportunities given all same cell opportunities are used. 
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Since only one of these three opportunities can used is a given 
cell in a timeslot, the maximum total sum is fixed and hence 

ci +C2 + C3 = C{p + q'). 

Define f{x)^EaM (ifr + which is simply the 

right hand side of ( 18 1. Because e > f{x), and because x G 
Hq n fli O ft2 n il^, we have: 



e > inf f{x) 

a;er2onf2in02nf23 



(19) 



Furthermore, for any function g{x) such that g{x) < f{x) for 
all X, and for any set 17 that contains the set ^loC^VllC^Vl2C^^^„ 
we have: 



e > inf g{x) 



(20) 



This follows because the function to be minimized is smaller, 
and the infimum is taken over a less restrictive set. We now 
define four new constraint sets rio, l^i, ^^2, ^^3 as follows: 



xw_ 
Ri 



Ri 



XaO < Cl 



C2 



a>2 



3=S X 



Xio 

Ri 



Ri 



^XaO 



Xqi 
R2 



< Ci+ C2+C3 



It can be seen that each of rig, i^i, 1^2! f^a is a subset of 
ClQ,Cti,Ct2,Ct2. Therefore, fig n fii H $^2 H il^ is a subset 



of fig n fii n 1^2 n rJa. Note that since 
following: 

12 12 

— < — < — < — 
Ri Ri i?2 R2 



< 



we have the 



(21) 



We now compute four different bounds for e, each having the 
form e > aX + (3. These bounds define the four piecewise 
linear regions of $(A). 

1) First note that f{x) > ^J2a,b^- This follows 
from the definition of f{x). Therefore taking g{x) = 
^a.b if' we have: 

_ . . r i sr^ Xab 

e > ml — > — — 

Because fig is given by b ^ab — NX, the above 
infimum is equal to Thus, we have our first linear 
constraint for any algorithm that yields a time average 
energy of e: 

e > ^ (22) 



Ri 



2) Next note that f{x) > 



xia 
NB.1 



This 



2 Xab 

Rl a,b ]y ■ 

ji^ T ^ > ^ for any non-negative integer 
pair (a, h) such that (a, &) V {(0, 0), (1, 0)} (using 



is because « 1 b ^ 2 



Therefore, taking this lower bound of f{x) as g{x), we 
have: 



e > inf 

cceOonOi 



JElO_ 



2 

111 



E 

(q,6)#(1,0) 



17 



The right hand side is equal to the solution of the 



following: 

Minimize: 

Subject to: 



xw 



2 

R~i 



E 

a,b 
(a,6)#(l,0) 



•^ab 



a,b 

Ri 



Xab = NX 



< Cl 



The above optimization is equivalent to minimizing 

WrI + mh:(^^~^w) subject to ^ < ci. The solution 



is clearly to choose Xio 
2A Cl 



e > 



i?iCi, and hence we have 
q 2 / 2Riq 



N 



R 



( ■ZHiq\ 
iT^^J 



(23) 



3) Next we have 

/(^)> 



Xw 

NRi 



2 \ ^ XaO 

R~i^W 

^ a>2 



1 \ - Xab 

R2^ N 



,6 
b#0 



which follows from the definition of f{x) and because 
^ < ^ for all positive 6 > 1. Thus, taking this lower 
bound of f{x) as g{x), we have: 



e > inf 

2:gf2onr2inn2 



JCiO_ 

NRi 



2 XaQ 

~Ri^l^ 

^ a>2 



R 



-E 



a,b 



•^ab 



This is equivalent to the following minimization: 



Minimize: 



Xw 

NRi 



NRi ^ 

^ a>2 



XaO 



1 



Subject to: 



NR2 

Xw 

Ri 



NX ~ Xw 



< Cl 



a>2 



XaO 



Xw 

Ri 



Ri 



XaO < Cl + C2 



a>2 



where we have aggregated the constraint J2a b ^ab = NX 
into the objective. The coefficients multiplying x^q and 
Sa>2 ■''ao are both negative, so that the above optimiza- 
tion is solved when xig + '2J2a>2^aO = Ri{ci + C2). 
Similarly, it can be shown that above optimization is 
solved when xw — Rici- This yields: 

A ^ (ci + C2) _ Ri 
R2 
P 



e > 



1 

9 ' Ra 
4) Finally, note that 

Xw ^ 
NRi Ri 



N 
X - 



NR2 

Rijp + q) 

29 



Cl 



C2 

2 



/(^)> 



E 



XaO 

N 



xoi 
NR2 



2 



E 

6>2 



(24) 



■^ab 



which follows from the definition of f{x) as well as 
for all b > 2. Taking this lower bound 



because 



< 



R2 
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of f{x) as g{x), we have: 

XlQ 



e > inf 



NRi Ri^ N 

^ ^ a>2 



E 

b>2 



where (l — (Iq O (li D CI2 D Cl^. This is equivalent to the 
following minimization (using J2a b ^ab = NX): 



Minimize: 



JElO_ 

NR^ 



NRi 



NRo 



NR- 



■ (^NX - Xio - ^ XaO - Xqi 



Subject to: 



Xio 

Ri 

XlQ 



< Cl 

2 



Rl i?i 



a;io 



^ XaO < Cl + C2 



a>2 



E^"0 



xai 

ID ■ D / , -aO -I- -5- < Cl + C2 + C3 

Ki Hi ^ n2 

a>2 



Letting y = J2a>2 and simplifying the optimization 
metric, the above optimization is equivalent to: 



Minimize: 



Subject to: 



XlO 

N 



Vi?i R2J 



y_(l 



xm 
NRo 



2A 

R^ 



XlQ 



< Cl 



Rl 

a^Lo _^ 2y 

Rl Rl 

a^LO ^ 2y ^ a^i 

Rl Rl i?2 
The above optimization is solved when xiq — RiCi, 
xin + 2y = Ri{ci + C2) and xoi = i?2C3. We thus have: 



< Cl + C2 

< Cl + C2 + C3 



e > 



2A , (C1+C2) 



R2 

p + q' 



N 
2 

R~2 



A- 



NR2^^''+''^-N 
Rljp + q) + 2R2q' 
29 



(25) 



The necessary set of conditions for $(A) function are obtained 
by combining these four bounds. ■ 

B. Proof of Sufficiency 

Now we present an algorithm that makes stationary, ran- 
domized scheduling decisions independent of the actual queue 
backlog values and show that for any feasible input rate A < /i, 
its average energy cost can be pushed arbitrarily close to the 
minimum value $(A) with bounded delay. However, the delay 
bound grows asymptotically as the average energy is pushed 
closer to the minimum value. Similar to the capacity achieving 
2-Hop Relay Algorithm, this algorithm also restricts packets to 
at most 2 hops. However, the difference lies in that it greedily 
chooses transmission opportunities involving smaller energy 
cost over other higher cost opportunities. An opportunity with 
higher cost is used only when the given input rate cannot 
be supported using all of the low cost opportunities. Thus, 
depending on the input rate A, the algorithm uses only a subset 



of the transmission opportunities as follows. 

1) If < A < all packets are sent using only source- 
destination transmission opportunities in the same cell. 

2) If ^ < A < all packets ai-e sent either using 
source-destination transmission opportunities in the same 
cell or source-relay and relay-destination transmission 
opportunities in the same cell. 

3) If ^^i^ < A < ^■i(p+gH2fl2g' ^ ^jj p^^]^gtg ggjjj 

using same cell transmissions (in either direct transmis- 
sion or relay modes), or adjacent cell source-destination 
transmission opportunities. 

4) And finally, when ^'^(p+i)+2^'^i' < A < /i, all transmis- 
sion opportunities that restrict packets to at most 2 hops 
are used. 

To make the presentation simpler, in the following, we only 
discuss the case where < A < ^^^-^j^-^- The basic idea 
and performance analysis for the other cases are similar. 

Let A = ^ + p^^^^^p^ where < p < 1 is a given 
constant. Also define a control parameter f3 (where 1 < /3 < 
1/p) that is input to the algorithm. This parameter affects an 
energy-delay tradeoff as shown in Theorem]?] 

Minimum Energy Algorithm: Every timeslot, for all cells, 
do the following: 

1) If there exists a source-destination pair in the cell, ran- 
domly choose such a pair (uniformly over all such pairs in 
the cell). If the source has new packets for the destination, 
transmit at rate Ri. Else remain idle. 

2) If there is no source-destination pair in the cell but there 
are at least 2 users in the cell, then with probability 
decide to use the same cell relay transmission opportunity 
as described in the next step. Else remain idle. 

3) If decide to use the same cell relay transmission oppor- 
tunity in step (2), randomly designate one user as the 
sender and another as the receiver. Then with probability 

(where < S < 1 and 6 — S{(3) to be determined 
later) perform the first action below. Else, perform the 
second. 

a) Send new Relay packets in same cell: If the transmitter 
has new packets for its destination, transmit at rate Ri. 
Else remain idle. 

b) Send Relay packets to their Destination in same cell: 
If the transmitter has packets for the receiver, transmit 
at rate Ri. Else remain idle. 

Note that the above algorithm does not use any adjacent cell 
transmission opportunities. All packets are sent over at most 
2 hops using only same cell transmissions. We now analyze 
the performance of this algorithm. 

Theorem 4: For the cell partitioned network (with N nodes 
and C cells) as described in Sec. ]Il] with minimum energy 
function $(A) as described above, and user mobility model 
as described in Sec. |II-B[ the average energy cost e of the 
Minimum Energy Algorithm with input rate A for each user 



such that A 



Rig 



P 



29 



(where < p < 1), a control 



parameter (3 (where 1 < (3 < 1 / p), and with S = 



$(A) + (;3-l)p( 



p-q 



satisfies: 



(26) 
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and the average packet delay D satisfies: 

ABNe{2d+l) 



D < 



(27) 



where i? is a constant given by (111 and d is a finite integer 
that is related to the mixing time of the joint user mobility 

, ( 4W^(p + i;)o/3 
r " V (P-<j)p(/9-l) 

process and is given by d = \og{i/i) — 

From the above, it can be seen that the control parameter 
/3 allows a {0{p - l),0(l/(/3 - 1))) tradeoff between the 
average energy cost and the average delay bound. Specifically, 
the average energy cost e can be pushed arbitrarily close to 
$(A) by pushing (5 closer to 1. However, this increases the 
bound on ~D as 1/{I3 - 1). 

Proof: The proof is similar to the proof of Theorem |2] 
and is given in Appendix C. ■ 

V. Capacity Gains by Network Coding 

Here, we show an example where the network capacity can 
be strictly improved by making use of network coding in con- 
junction with the wireless broadcast advantage. Specifically, 
consider a network with 6 nodes and 4 cells. Suppose the 
steady-state location distribution for all nodes is uniform over 
all cells. Thus, tt^ = 1/4 for all c. The one-to-one traffic 
pairing is given by 1 o 2, 3 -H> 4, 5 -;-?> 6. Let Ri = 1 
and i?2 = 0. Thus, this example only allows same cell 
transmissions. We further assume that when a node in a cell 
transmits, all other nodes in that cell receive that packet. Note 



that the 2-Hop Relay Algorithm presented in Sec. III-B does 
not make use of this feature. 

Using Theorem T] the network capacity under the model 
presented in Sec. |ll can be computed. Specifically, the network 
capacity is given by /i = packets/slot per node where 



and using (2i, we have q = 1 — {\ 



and p 



We now show how network coding can be used to achieve 
a throughput that is strictly higher than /i. First we define 4 
distinct configurations of the nodes. In configuration I, nodes 
1,4, and 5 are in the same cell and the other nodes can be 
in any of the remaining cells (but not in the same cell as 
nodes 1,4, and 5). Note that this cell can be any one of 
the 4 cells. From the assumption about the node mobility 
process, the steady-state probability of configuration I is given 
by i^44 X [\Y X (|)^. In configuration II, nodes 2,3, and 5 
are in the same cell and the other nodes can be in any of the 
remaining cells (but not in the same cell as nodes 2. 3, and 5). 
In configuration III, nodes 2, 4, and 5 are in the same cell and 
the other nodes can be in any of the remaining cells (but not 
in the same cell as nodes 2, 4, and 5). Finally, in configuration 
IV, nodes 1, 3, and 5 are in the same cell and the other nodes 
could be in any of the remaining cells (but not in the same 
cell as nodes 1, 3, and 5). Note that these configurations cannot 
occur simultaneously as each consists of node 5. Further, the 
steady-state probability of each configuration is given by v. 

In the following, we will modify the 2-Hop Relay Algorithm 



and demonstrate an improvement in the throughput of nodes 
1, 2, 3 and 4 over /i. For each configuration, we will only focus 
on the transmissions in the cell with the three nodes that define 
that configuration. The 2-Hop Relay Algorithm for the other 
cells remains the same. 

Note that under each configuration, there are no source- 
destination pairs in the cell of interest. Thus, under the 2-Hop 
Relay Algorithm, a node is selected as the transmitter with 
probability | while the remaining two nodes are equally likely 
to be selected as the receiver. Further, the transmitter is sched- 
uled to transmit a new packet to the receiver with probability 
and is scheduled to transmit a relay packet to the receiver 
with probability Thus, in each configuration, each of the 
two nodes other than node 5 is scheduled to transmit a new 
packet to node 5 with probability 5 x ^ x ^^^^^ = ^^^pp-- 
Also, in each configuration, node 5 is scheduled to transmit a 
relay packet to each of the other two nodes in the cell with 



probability i x i x 



Adding the probabilities 



3 2^2 12 

associated with these four scheduling decisions yields 



(l-S) , (l-S) , (l + S) , (1 + 5) 



12 



12 



12 



12 



(28) 



We now modify the 2-Hop Relay Algorithm to take ad- 
vantage of network coding. For all configurations other than 
the four as defined above, the algorithm remains the same. 
However, in each of the configurations I, II, III, IV, we change 
the probability of scheduling a node to transmit a new packet 
(for relaying) to node 5 from ^^jp to | x -^3-^ — 
where < e < 1. Also, node 5 is scheduled to transmit a 
relay packet to the other two nodes in the cell with probability 
i X ^^Ij^'^-' — However, whenever node 5 has at least 

one packet for each of the two other nodes, it broadcasts 
a XOR of two packets destined for these nodes in a single 
transmission. If node 5 does not have at least one packet for 
each of the two other nodes, it would simply transmit a regular 
packet (if available). Note that under the original 2-Hop Relay 
Algorithm, the two scheduling decisions of node 5 transmitting 
a relay packet to the other two nodes are taken with probability 
^"^jj^ each. These are now replaced by a single scheduling 
decision of node 5 broadcasting a XORed relay packet and 
this has probability The probabilities associated with 

the other scheduling decisions under this modified algorithm 
remain the same as the original 2-Hop Relay Algorithm. The 
sum of probabilities associated with the modified scheduling 
decisions as described above is given by 



(1-e) ^ (1-e) , (l + 2e) 



(29) 



of Sec. III-B when one of these configurations occur in any cell 



9 9 9 3 

This is the same as ( |28] l. Thus, it can be seen that the 
probabilities of all scheduling decisions under the modified 
algorithm sum to 1. 

To see how the nodes can recover the original packets 
from the XORed packet, we further classify each configuration 
into type A, B and C depending on the scheduling decision 
as shown in Fig. [3] The configurations of type A and B 
correspond to the scheduling decisions in which a node is 
scheduled to transmit a new packet (for relaying) to node 5. 
The configurations of type C correspond to the scheduling 
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Fig. 3. An example showing capacity gains possible by using network coding Fig. 4. Additional relay queues at node 5 under the network coding enhanced 
in conjunction with the wireless broadcast advantage. 2-Hop Relay Algorithm that is used in configurations I, II, III, IV. 



decisions in which node 5 is scheduled to transmit relay 
packets to the other two nodes (either as a network coded 
XOR packet whenever possible or a regular packet). In each 
configuration of type A or B, whenever a new packet is 
transmitted by a node to node 5 for relaying, the other node 
overhears the packet and stores a copy. For example, in Fig. 
I-A, when node 1 transmits a new packet (destined for node 
2) to node 5, node 4 overhears this transmission and stores a 
copy of this packet. Similarly, in II- A, when node 3 transmits a 
new packet (destined for node 4) to node 5, node 2 overhears 
this transmission and stores a copy of this packet. In each 
configuration of type C, whenever node 5 has at least one 
packet for each of the two other nodes, note that each of 
these two nodes already has a copy of the packet destined 
for the other node (that it obtained by overhearing earlier 
in a type A or B configuration). Therefore, when node 5 
transmits a XOR packet, both of these nodes can recover 
the original packets destined for them by using the side 
information already available to them in the form of previously 
overheard and stored packets. For example, in III-C, node 5 
is in the same cell as nodes 2 and 4 and suppose it has a 
packet for each of them. Then, node 2 must have the packet 
that is destined for node 4 that it overheard in II-A. Similarly, 
node 4 must have the packet that is destined for node 2 that 
it overheard in I-A. Thus, when node 5 broadcasts a XOR 
packet in a single transmission, both nodes 2 and 4 can retrieve 
their desired packets. Thus, this single transmission effectively 
delivers two packets. Note that under a scheme that does not 
allow mixing of packets, at most one packet can be transmitted 
per transmission. 

To demonstrate gains in throughput under this "network 
coding enhanced" 2-Hop Relay Algorithm, we define the 
following additional relay queues at node 5 as shown in Fig. 
|4] Arrivals to and departures from these queues happen only 
when scheduling decisions corresponding to the 12 configu- 



rations in Fig. [3] are made according to the enhanced 2-Hop 



Relay Algorithm. U-j\t) and U-j'{t) refer to the queue of 
packets destined for nodes i and j respectively that will be net- 
work coded whenever possible. Fig. |4] shows the arrival rates 
and the corresponding configurations (when arrivals happen to 
these queues) as well as the service rates and corresponding 
configurations (when packets are served from these queues). 
Note that each queue has an arrival rate of ^^^g^"^ and sees 
a service rate of Since (1 + 2e) > (1 — e), all these 

queues are stable. The additional throughput for nodes 1,2,3 
and 4 over the 2-Hop Relay Algorithm without network coding 



is given by 



2(l-e) 2(l-(5) 



9 



1/ packets/slot. This is strictly 

For example, by choosing 
packets/slot is achievable. 



36 



12 

positive for any < e < 
e = |, a throughput gain of 
Thus, the capacity can be strictly increased over a scheme 
that is restricted to pure routing. 

VI. Conclusions 

In this work, we investigated two quantities of fundamen- 
tal interest in a delay-tolerant mobile ad hoc network: the 
network capacity and the minimum energy function. Using 
a cell-partitioned model of the network, we obtained exact 
expressions for both these quantities in terms of the network 
parameters (number of nodes and number of cells C) and 
the steady-state location distribution of the mobility process. 
Our results hold for general mobility processes (possibly non- 
uniform and non-i.i.d.) and our analytical technique can be ex- 
tended to other models with additional scheduling constraints. 

We also proposed two simple scheduling strategies that 
can achieve these bounds arbitrarily closely at the cost of 
an increased delay. Both these schemes restrict packets to 
at most 2 hops and make scheduling decisions purely based 
on the current user locations and independent of the actual 
queue backlogs. For both schemes, we computed bounds on 
the average packet delay using a Lyapunov drift technique. 
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In this paper, we have focused on network control al- 
gorithms that operate according to the network structure as 
presented in Sec. We assumed that the packets themselves 
are kept intact and are not combined or network coded. As 
shown in the example in Sec. [V] it is possible to increase the 
network capacity by making use of network coding and the 
wireless broadcast feature. An interesting future direction of 
this research is to determine the exact capacity region with 
such enhanced control options. 

Appendix A 
Proof of Lemma 1 

Here, we prove the bound in Lemma [T] We have 
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where step two follows from the independence of node mo- 
bility processes, step three follows from ([T]i and the second 
last step uses the inequality (1 — aj'^)'^ > (1 ~ 2Na'^'^). 
This can be shown by induction as follows. This holds for 
iV = 1. Suppose it holds for some integer i > I, i.e., 
(1 — a^'^y > (1 — 2ia'f'^). Then, we have 

(1 - aj'^y+^ = (1 - a7'^)''(l - aj'^) > (1 - 2iaj'^){l - aj"^] 

> (1 - 2(^ + l)a7'') 

The upper bound can be shown similarly, except that we use 
the inequality (1 + 07'*)^ < (1 + 2Naj'^) for all iV > 2 
whenever d is such that aj'^ < To show this, let a-f"^ = 

c/N^ where < c < 1. Note that < c/N < 1. Then 

N{N -1) 



in the steady-state. Then the expected number of cells with 
a source-destination pair is given by E |^^-^ By 

linearity of expectations, this is equal to X^ljLi E {/c(^)}- 
To compute E{/c(i)} for any cell c, note that by the inde- 
pendence of user mobility processes, tt^ is the probability 
of finding any particular source-destination pair in cell c in 
the steady-state. Since there are N/2 such pairs and they 
occur independent of each other, the probabiUty of finding 
no source-destination pair in cell c in the steady-state is given 
by (1 — TT^)^/^. Thus, the probability of finding at least 1 
source-destination pair is 1- (l-7r2)^/2. Using this, we get 

'z = ^Ef=i(i-(i--cW. 

Derivation of p: To compute the probability of finding at 
least 2 users in a cell c, we note that this can be obtained by 
first computing the probabilities of finding no user and exactly 
1 user in cell c and then subtracting these from 1. These are 
given by (1 — tt^)^ and (^)7rc(l — tTc)^^^ respectively. Using 

this, we get p - i Ec^ii^ - (1 - ^c)^ - NtTc{1 - ^c)^"^). 

Derivation of q' : The probability of finding exactly 1 user in 
cell c is given by (^)7rc(l-7rc)^~^ The probability of finding 
its destination in an adjacent cell given that it is not it cell c is 
given by Tli^B '^i which we have defined as Jladj{c)- 
Using this, we get 9' = ^ Ef=i(nadj(c)iV^c(l - ^c)"^"')- 

Derivation of p' : Given that there is exactly 1 user in cell 
c, the probability that at least 1 of the remaining iV — 1 users 
is in an adjacent cell is given by 1 — (1 — liadjic))^^^ ■ Thus, 
wegety = iEti(l-(l-nad,(c))~'i)iV7r,(l-^,)^-i. 

Derivation ofq": We first compute the probability of finding 
i users in cell c such that there are no source-destination pairs. 
Clearly, 1 < i < ^ since there must be at least 1 source- 

(N/2\ 

destination pair for i > ^. Next, note that 2' is the 

(J 

probability of finding no source-destination pair in a cell given 
that there are i users in that cell. (^)7r*(l — tTc)^^' is the 
probability of having i users in cell c. Finally, the probability 
that there is at least 1 node in an adjacent cell that will make 
a source-destination pair with one of these i nodes given that 
it is not in cell c is given by (1 — (1 — ^adj{c)y). Combining 
all these, we get 

C N/2 



= 1 

< 1 

< 1 



d\2 



C 

N 



N - 1 



= 1 



N 



< 



1 N 



1 - 



N 



N-l - ' 



N-c 
yN>2 



Appendix B 
Derivation of Probability Expressions 

In what follows, we will use the linearity of expectations 
property to compute the probability expressions in Q. 

Derivation of q: Let Ic{t) be an indicator variable that 
is 1 if there is a source-destination pair in cell c in slot t 



c=l i=l ^ ^ 

Derivation of p" : Similar to the derivation of q", the 
probability of finding i users in cell c such that there are no 
source-destination pairs in cell c as well as any adjacent cells 
is given by 2^(^/2)7r^(l - - Iladjic)y. Since we 

also want at least 2 users in cell c, we sum from i = 2 to 
N/2. This yields 

c=l i=2 V * / 



Appendix C 
Proof of Theorem|4] 

Here, we establish the bounds (|26]l and dZT] 
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When ^ < A < ^^^^j^, under the Minimum Energy 
Algorithm, all transmissions are either same cell direct trans- 
missions or same cell relay transmissions. Specifically, each 
user either transmits directly to its destination or transmits 
new packets to a relay or transmits relayed packets to their 
destination in the same cell. Each such transmission involves 
one unit of energy cost and therefore the average energy 
cost per user e can be expressed in terms of the rates of 
these transmission opportunities. The rate at which same cell 
direct transmissions are scheduled is given by Cq. The rate at 
which same cell relay transmissions are scheduled is given by 
PpC{p — q). Thus, we have: 



Cq 



PpCip^q) 



N 

$(A) + (/?-l)p 



Thus, e can be pushed arbitrarily close to $(A) by choosing 
/3 close to 1. 



The delay of the Minimum Energy Algorithm can be 
analyzed using a procedure similar to the one used in the 
proof of Theorem|2] We first evaluate bounds on the expression 
in ( [12I by computing the steady-state service rates Jl]^ jjl^^} 
achieved by the Minimum Energy Algorithm. We have the 
following 2 cases: 



1) Node i generates type c packets: In this case, 
^[A'f\t)] = A and E./^if (*) = 0. To calculate E^A^L'^ 
define ri,r2,r3 similar to that in the proof of Theorem 
|2] Then, the total rate of transmission over the network 
is given by N{ri + r2 + r^). Similar to Theorem |2] we 
have r2 = ytI^'s- Since only same cell direct transmissions 
are used, we have Nri = CRiq. Also, a same cell relay 
transmission is scheduled with probability /3p whenever there 
is no source-destination pair in the cell but there are at least 
2 users in the cell. Thus, the sum total transmissions in the 
network can be expressed in terms of the quantities p and q 
as N{ri +r2+ rg) = C{Riq + i?i/3p(p - q)). Solving for 
ri,r2,r3, we have: 



Riq Ri{p-q){l-5)Pp 



261 



Ri{p-q){l + S)pp 



Therefore, we have: 



E 



-(c) 

P'ib 



ri + r2 



29 



Riq , Ri{p~q){l-5)(ip 



(30) 



Let 5 = ^ and aj'^ 



{p-q)pS 



29 



{p-q)p{fi-l) 



Note that this choice of S can be shown to represent a valid 
probability, because 1<P<^=>0<^ - ^ < ^ - f =^ 



< (5 < 1. Then, using (|T3]l, the last term of (12i under this 



case can be expressed as: 



(ri + r2)(l - 2Na-f'^) - X > {n + ra) 



20 



(1-^)^-^-1 



> 



_ Ri{p~q)p5 
29N 

R,{p-q)p{l3-l) 



A 



^1 + t™?' (^1 + 



where we used the relations A — ; 1 r (jTrjyg- 
r2)2NaY < ^^^^^ and (1 - <5)/3 - | - 1 > 
These can be shown as follows. Using pO| , we have A — 

Rii ^ ^ Rijp-i) - ^, J r:2_ p^gxt; 



P- 



28 



ri 



(l-5)/3- 



^ Ri{p~q)p5 
29N 

Finally, using 5 = 4?^, we have: 



213 



/3+1 



N 



5 
N 



1 > 



13-1 



2 4^-4 



2) Node i relays type c packets: From our traffic model, 
we know that in this case A['^\t) = for all t. To compute 

l^b'Pib' '^^^ 'l^a'PaA^ ^'^^^ ^^^^ Minimum Energy Algo- 
rithm schedules relay transmissions such that all — 2 relay 
packet types are equally likely. Thus we have: 



b 



iV- 2' 



E 



-(c) 
Pai 



r2 



Let 6 4r7r and 017'^ 



{p-q)pS 



N 



{p-q)p(f}-l) 



213 ~ 2{p+q)N^ ~ 48(v+a)N^ ^ ■ 

Then, using ( pj) , ( [l4| , the last term of ([12]) under this case 
can be expressed as: 

> ( E (1 - 2 - ( E (1 + 2^"^') 



> 



^3 - r2 
N-2 
Rl{p-q)pP 
9N 



ra + r2 
N -2 J 



\ 2Na-/'^ 



N 



> 



fli(p-g)p(/3-l) 
4iV6l 



where we used the inequality 2Na^ < Combining 
these two cases, with 6 = and 07^^ — ''^^l^p'^^j^^^V 



we 

Rlip-q)p(j3-l) 



have E {e, pl'.^t) - t^ilit) - A^\t)\Uit - d)} 
~4Ne ^~^^ • Using this in |l2| , we get: 
A{t,d) <BN{2d+l) 

i—l c^i 



> 



4A^6l 
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This is in a form that fits (|6]l. Using the Lyapunov Drift 
Lemma, we get 



limsup-> > E\U}'(t)}< — ^— — 

i^ofc ^ ^ - Riip-q)pW-l) 

The total input rate into the network is A^A. Thus, using 
Little's Theorem, the average delay per packet is bounded by 

4:BNe{2d+l) 
\Ri{p-q}p{fl-l}- 
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